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DYNAMIC RATE CONTROL METHODS AND APPARATUS FOR SCHEDULING 
DATA TRANSMISSIONS IN A COMMUNICATION NETWORK 

5 Related Applications 

The present invention is related to the inventions described in U.S. Patent AppUcation Serial 
No. 09/517,659 filed in the name of inventors S. Stolyar et al and entitled "Method of Packet 
Scheduling, With Improved Delay Performance, for Wireless Networks," and U.S. Patent 
AppUcation Serial No. 09/393,949 filed September 10, 1999 in the name of inventors K. Ramanan 
10 et al. and entitled "Method and Apparatus for Scheduling Traffic to Meet Quahty of Service 
; -| Requirements in a Communication Network," both of which are commonly assigned herewith and 
15; incorporated by reference herein. 

m Field of the Invention 

The present invention relates generally to wireless communication systems and more 
particularly to techniques for scheduling data transmissions in wireless code division multiple access 
Q (CDMA) wireless systems and other types of communication networks. 

Q Background of the Invention 

20 Next-generation wireless networks are expected to support a wide range of services, 

including high-rate data apphcations. Examples of such wireless networks include those CDMA 
systems known in the art as IS-95B, cdma2000, third generation (3G), Universal Mobile Telephone 
System (UMTS) and wideband CDMA (WCDMA) systems. The high-rate data applications to be 
supported by these networks are generally characterized by bursty traffic. Moreover, in contrast to 

25 voice users, data applications can usually accept a certain amount of packet delay, as long as 
sufficient data throughput is provided over a particular interval. As a result of the bursty traffic and 
relative delay intolerance of data apphcations, scheduling techniques can be particularly effective 
in optimizing system throughput for such applications. 

An example of a conventional scheduling technique particularly directed to the scheduling 

30 of high-rate data transmissions on a base-to-mobile hnk, i.e., downlink, of a CDMA system is 
described in A. Bedekar et al., "Downhnk scheduling in CDMA data networks," Proc. IEEE 
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Globecom '99, pp. 2653-2657, 1 999, which is incorporated by reference herein. Related techniques 
may be applied in low-mobility scenarios, such as indoor networks. In such environments, Rayleigh 
fading frequencies can be quite low, and the fading levels can even be anticipated to some extent. 
For example, fading can be measured by having the base station provide a pilot signal which can be 
measured by all the users and utihzed to generate corresponding signal-to-noise ratios (SNRs). The 
measurements can be fed back to the base station, and used to estimate fading levels and SNRs, and 
hence user rates, in subsequent slots. It is therefore possible to exploit the variations in the channel 
conditions so as to transmit to the user with the currently "best" channel. This type of an approach 
is utilized in the High Data Rate (HDR) technique developed by Qualcomm, More particularly, in 
such a technique mobile users may be configured to make measurements on a common pilot signal, 
and to report the measurements on a dedicated uplink channel, e.g., a 3.2 kbps dedicated uplink 
channel. These and other aspects of the HDR technique are described in greater detail in P. Bender 
et al., "CDMA/HDR: a bandwidth-efficient high-speed wireless data service for nomadic users," 
IEEE Communications Magazine, Vol. 38, pp. 70-77, 2000, which is incorporated by reference 
herein. 

The above-noted estimates of user rates for subsequent slots can provide a basis for 
scheduling user transmissions. However, the estimates fail to provide the probability that each user 
will have a given SNR and a potential maximum rate. Although it is possible to determine such 
probabilities from a joint distribution of user rate declarations, this would generally involve keeping 
records from slot to slot of the user rate declarations, and building up a relative frequency table, and 
would thus require excessive computational and memory resources in a practical system 
implementation. A need therefore exists for improved techniques for scheduling high-rate data 
transmissions, on a CDMA downlink or other type of communication network link, without 
requiring direct estimation of the frequency of occurrence of particular user rates. 

Summary of the Invention 

The invention provides methods and apparatus for scheduling data transmissions in CDMA 
systems and other types of communication networks. 
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In accordance with one aspect of the invention, the scheduhng of data transmissions for a 
CDMA system downhnk or other type of communication network is implemented using a revenue- 
based policy. For a given transmission slot or other transmission interval, a maximum-rate user is 
identified from among a set of users requesting data transmissions, and a data transmission of the 
maximum-rate user is scheduled for the given interval. The maximum-rate user is identified based 
on apphcation of coefficients of a revenue vector to corresponding feasible rates of the requesting 
users. The revenue vector is determined in an iterative manner using an adaptive algorithm which 
updates the revenue vector periodically to compensate for observed deviations between actual and 
target throughput. As a result, the deviations are reduced over time and the revenue vector converges 
to an optimal revenue vector. Advantageously, the revenue vector can be determined without direct 
estimation of the frequency of occurrence of particular user rates. 

In accordance with another aspect of the invention, the maximum-rate user for which a data 
transmission is scheduled in an n-th transmission interval is a user m\n) identified as follows: 

m'^in)^ arg max w^R^ (n) , 

where Wj, . . . , denote the coefficients of a revenue vector w, and RJn) denotes the feasible rate 
for an m-Xh. one of M users for the n-th transmission interval. 

In accordance with a further aspect of the invention, convergence of the revenue vector to 
its optimal value is ensured by the adaptive algorithm periodically adjusting the coefficients of the 
revenue vector in accordance with a specified step size taken from a predetermined convergent 
sequence. 

In accordance with yet another aspect of the invention, the adaptive algorithm updates the 
revenue vector based on sample periods which increase in size as a function of time. For example, 
the updates to the revenue vector occur at predetermined slots K{n), with the length of an n-th sample 
period given by L{n) := K(n + 1) - K{n), The adaptive algorithm may be configured to update the 
revenue vector w(n) for an n-th sample period in accordance with an update direction v(w(n)). The 
update direction is determined based on throughput obtained during the n-th sample period when the 



Borst 13-7 

revenue vector w(n) is used. At a subsequent update of the revenue vector for an (n + l)-th sample 
period, the revenue vector is recursively determined as 

w(n + 1) = w(n) - ?>{n)v{w{n)), 

where 5(n) denotes a step size of the above-noted predetermined convergent sequence. The update 
direction v{w{n)) may be determined by averaging the actual throughput for each user over the 
sample period, partitioning the users into a first group with above-average throughput and a second 
group with below-average throughput, decreasing the revenue vector coefficient values for the users 
in the first group, and increasing the revenue vector coefficient values for the users in the second 
group. As the size of the sample period increases, the actual throughputs will approach specified 
target throughputs. 

Brief Description of the Drawings 

FIG. 1 shows an exemplary code division multiple access (CDMA) wireless communication 
system in which a technique for scheduling of data transmissions in accordance with the invention 
is implemented. 

FIG. 2 is a block diagram showing a base station which is utilized in the FIG. 1 system for 
scheduling data transmissions in accordance with the invention. 

FIG. 3 is a flow diagram of the a data transmission scheduling process that is implemented 
in the base station of FIG. 2. 

FIG. 4 shows a plot of expected throughput per time slot as a fimction of revenue vector 
coefficient for an example two-user case. 

Detailed Description of the Invention 

The present invention will be illustrated below in conjimction with an exemplary code 
division multiple access (CDMA) wireless communication system. It should be understood, 
however, that the invention is not limited to use with any particular type of communication system, 
but is instead more generally appUcable to any wireless system in which it is desirable to provide 
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improved scheduling of data transmissions. For example, although the techniques are illustrated 
with reference to certain types of CDMA systems, it will be apparent to those skilled in the art that 
the techniques are also applicable to other CDMA systems, as well as to other types of wireless 
systems and communication networks. 

Example System 

FIG. 1 shows an exemplary CDMA wireless cellular system 1 0 in which a data transmission 
scheduling technique of the invention is implemented. The system 10 may be configured in 
accordance with well-known conventional wireless communication system standards. The system 
10 as illustrated includes a mobile station (MS) 12 and a number of base stations BSl, BS2, BS3 and 
BS4. The base station BSl of FIG. 1 represents a primary base station, communicating with mobile 
station 12 via the path designated by solid two-way arrow 24, while the surrounding base stations 
BS2, BS3 and BS4 may serve as secondary base stations which can detect signals transmitted via 
the paths indicated by dashed one-way arrows 26. The primary base station BS 1 communicates with 
the mobile station 12 using CDMA techniques that are well known in the art. As the mobile station 
12 moves within the system 10, handoffs occur, such that base stations other than BSl become 
primary base stations for communicating with the mobile station 12. Also, handoffs may occur from 
a given antenna sector of the primary base station and another sector of that base station. 

Transmissions from a given base station of the system 10 to the mobile station 12 are referred 
to as base- to-mobile or downlink transmissions. Transmissions from the mobile station 12 to a given 
base station of the system 1 0 are referred to as mobile-to-base or uplink transmissions. 

The system 10 in this illustrative embodiment also includes first and second mobile switching 
centers (MSCs) 14-1 and 14-2. A given MSC typically connects several BSs with apubhc switched 
telephone network (PSTN) 16. For example, MSC 14-1 connects base stations BSl and BS2 with 
the PSTN 16 and MSC 14-2 connects base stations BS3 and BS4 with the PSTN 16. The system 
1 0 also includes a memory 1 8 having a number of registers including a home location register (HLR) 
20 and a visitor location register (VLR) 22. The HLR 20 and VLR 22 store user data and billing 
information for each mobile station 12 of the system 10. 
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The system of FIG. 1 is shown in simphfied form for purposes of illustration, and will 
generally include a substantially larger number of base stations and mobile stations, as well as other 
system elements that are not explicitly shown in the figure. Each of the mobile stations within the 
system is more generally referred to herein as a "user." It should be noted that the term "user" is 
intended to include not only mobile stations, but also any other type of subscriber unit of a given 
system, including, e.g., fixed wireless subscriber units. 

FIG. 2 shows a more detailed view of an apparatus which represents a given one of the base 
stations BSl, BS2, BS3 or BS4 of the system 10. The apparatus includes an antenna 32 which 
receives signals fi^om and transmits signals to other elements of the system 10. A receive signal is 
directed by a diplexer filter 33 to an input of a receiver 34 which may implement conventional 
downconversion, demodulation, digital-to-analog conversion and other processing of the receive 
signal. A transmitter 35 performs complementary operations such as analog-to-digital conversion, 
modulation and upconversion to generate a transmit signal which is directed via diplexer 33 to the 
antenna 32 for transmission. 

A processor 36 is coupled to both the receiver 34 and the transmitter 35. The processor 36 
operates in conjunction with a memory 3 8 to control the communication functions of the base station 
BSz. For example, data or other information in messages which are part of a receive signal may be 
supplied to or otherwise processed by a processor 36 such that processor 36 can implement data 
transmission scheduling techniques to be described in greater detail below. 

General Scheduling Approach 

The present invention provides improved techniques for scheduhng of high-rate downlink 
data transmissions in a CDMA wireless system such as that illustrated in conjunction with FIGS. 
1 and 2. Advantageously, the techniques of the invention avoid the need for direct estimation of the 
fi-equency of occurrence of particular user rates. 

As previously noted, conventional scheduling techniques for high-rate data transmissions on 
a CDMA downhnk may be configured to estimate the channel conditions, and to transmit to the user 
with the currently "best" chaimel. In accordance with an illustrative embodiment of the invention, 
the user having the "best" channel may be identified as a maximum-rate user when a set of feasible 
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rates are weighed with appropriately-determined coefficients. Each of these coefficients may be 
viewed as a particular price or other reward value in a revenue vector, and the data transmission 
scheduling may then be viewed as a revenue-based pohcy which always assigns an available 
transmission slot to the user yielding the maximum revenue. The term "coefficient" as used herein 
in conjunction with a revenue vector is intended to include a price, reward value or any other 
revenue-related value that may comprise an element of the revenue vector. As will be described in 
detail below, the revenue vector coefficients are used to weight a set of feasible rates for system 
users ehgible for data transmission, and are adjusted over time such that actual throughput ratios 
converge to target throughput ratios. 

It should be noted that direct calculation of an optimal revenue vector of the type described 
above may be a difficult task, in that such direct calculation will generally require detailed 
information on the channel statistics. The present invention avoids such direct calculation by 
providing adaptive algorithms for determining the optimal revenue vector in an iterative manner, 
without the need for direct estimation of the frequency of occurrence of particular user rates. 
Starting firom an arbitrary initial vector, the algorithms iteratively adjust the revenue vector 
coefficients to compensate for observed deviations fi-om target throughput ratios. The adjustments 
ensure that discrepancies between actual and target throughput cannot persist, and are gradually 
reduced in size over time such that the revenue vector converges to the optimal revenue vector. 

The illustrative embodiment will be described in conjimction with a base station serving M 
data users. The base station transmits in slots of fixed duration in this embodiment, although 
transmission intervals other than fixed-duration slots may be used in other embodiments. In each 
slot, the base station transmits to exactly one of the users. It is assumed that the feasible rates for 
the various users vary over time according to a stationary discrete-time stochastic process 
. . . , Rm{^)}, h = 1, 2, . . withi2^(n) representing the feasible rate for user m in the n-\h slot. It is 
assumed that the base station has perfect knowledge of the maximum feasible rate (n) for user m 
at the start of the n-th slot. Let (i?^ , . . . , ) be a random vector with a distribution corresponding 
to the joint stationary distribution of the feasible rates. 

It should be noted that these and other assumptions made herein are for purposes of 
simphfication and clarification of the description, and should not be construed as limiting the scope 
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of the invention in any way. Alternative embodiments in which one or more of the assumptions do 
not apply will be apparent to those skilled in the art. 

The above description allows for dependence between the feasible rates for the various users. 
Independence may be a reasonable assumption in the case of an isolated base station serving a group 
of independent users. In the case of several base stations, however, the feasible rates may vary not 
only due to independent fading, but also because of the common impact of control actions at adjacent 
base stations. For example, base stations may transmit at reduced power if there are no backlogged 
users, inducing strong correlations in interference levels between users. 

The slot duration in the illustrative embodiment may be assumed to be relatively short 
compared to the relevant time scales in the traffic patterns and delay requirements of the data users. 
For example, the slot duration may be on the order of the 1.67 ms used in the above-noted HDR 
technique. 

The M data users in the illustrative embodiment may be viewed as a subset of active 
backlogged users among a greater population, which may change over time. For scheduling 
purposes, however, the separation of time scales allows the subset of active users to be viewed as 
nearly static and continuously backlogged. In practice, flow control algorithms such as 
Transmission Control Protocol (TCP) will typically be used to feed the data into the base station 
buffer at a relatively slow rate, comparable to the actual throughput provided to the user over the 
wireless link. Thus, the bulk of the backlogs will usually reside at the sender rather than the base 
station buffer. 

Throughput maximization is a common performance objective for high-rate data 
transmissions on a CDMA downlink, and may be the performance objective in the illustrative 
embodiment of the invention. It should be noted that this objective can be achieved simply by 
assigning each slot to the user with the currently highest feasible rate. The disadvantage of this 
approach is that typically only a few strong users will ever be selected for transmission, causing 
starvation of all others. An alternative approach which alleviates this problem is to equalize the 
expected throughput of the various users. This can easily be achieved by assigning each slot to the 
user with the currently smallest cumulative throughput. The disadvantage of this alternative 
approach is that it does not exploit the variations in the feasible rates. Moreover, by insisting on 

8 
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equal throughput, a few weak users may cause the throughput of all others to be dramatically 
reduced. Another approach is to equalize the proportion of slots allotted to the various users. Such 
an approach can be realized simply by using a round-robin scheduling technique. Again, however, 
this approach fails to take advantage of the fluctuations in the feasible rates. In addition, some users 
may end up with extremely low throughput, despite receiving their fair share of the number of 
transmission slots. 

The present invention in the illustrative embodiment avoids the problems of the above-noted 
approaches by defining throughput target values , . . . , for the various users. These target 
values may be set arbitrarily, taking into account the Quality-of-Service (QoS) requirements of the 
users, or possibly their current activity levels or locations. For example, the target values may be 
set lower for users with higher path losses, in order to prevent weak users from dragging down the 
throughput of all other users. The target values may also be applied to the proportion of transmission 
slots allotted to the various users. 

Some additional notation will now be introduced. Let Y^(n) :=XJ^n) RJn), withX^(n) a 0-1 
variable indicating whether or not the n-th slot is assigned to user m. Further define: 



as the expected average throughput received by user m after N slots. The illustrative embodiment 
considers the problem of maximizing the minimum long-run expected average normalized 
throughput 

niin V /oc^^ 



with := lim inf^^^j^^ {N), The above problem may equivalently be formulated as maximizing z 
subject to the constraint z < /a^ for all = 1, . . . , M . The constraints may in fact be tightened 
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to z^y^ /a^, since one can always reduce the throughput for the users with a surplus. With the 

M 

equaUty constraints in place, the obj ective function may then be generalized to ^ w^y^ / a ^ for any 

positive vector w gR^ . 

The above-stated problem is thus equivalent to maximizing an arbitrary weighted sum of the 
5 throughputs, subject to the normaUzed throughputs being equal. Thus, any scheduling strategy 
which maximizes an arbitrary weighted combination of the throughputs, while equalizing the 
normalized throughputs, is optimal. 

A class of revenue-based scheduling strategies for implementing the above-described 
maximization will now be introduced. Assume there are rewards w^, . . . , per bit transmitted 
i|0 to the various users. A revenue-based scheduUng strategy assigns the n-th transmission slot to the 
U| user rn(n) with the current maximum rate-reward product, i.e., 

m^{n)= argmax>i;^7?^(^). 

it|5 Clearly, the above revenue-based scheduling strategy maximizes the revenue eamed in each 

individual slot, and thus the total cumulative revenue as well as the average revenue. Different tie- 
breaking rules may be used. However, regardless of the particular tie breaking rule used in a given 
implementation, this revenue-based scheduling strategy will definitely not assign the n-th slot to any 
user k with 

20 

It should be noted that revenue is simply a weighted combination of the throughputs, and 
thus any revenue-based policy which balances the throughputs, is in fact optimal. It should also be 
25 noted that setting throughput targets is equivalent to normalizing the feasible rates by the 

10 
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corresponding values. Therefore, it will be assumed in the following description that the throughput 
target values are discounted for in the rates, such that (a^, . , . , a^) = (1, . . . , 1), 

In practice, there is always a small probabihty that a transmission fails because the signal 
cannot be successfully decoded. The results described herein remain valid if R^(n) is redefined to 
represent the expected feasible rate, and the 0 - 1 variable X^(?7) is modified to indicate both which 
user is selected and whether or not the transmission is successful 

FIG. 3 is a flow diagram illustrating the general data transmission scheduling process 
described above. This process is implemented in the base station of FIG. 2, e.g., using software 
stored in memory 38 and executed by processor 36. In step 100, the base station determines the 
feasible rates , . . . , i?^ ) for the M users eligible for data transmission in the n~th transmission 
slot. The base station then applies in step 102 the coefficients of the current revenue vector to the 
corresponding feasible rates, so as to generate a set of weighted feasible rates for the M users. In 
step 1 04, the base station identifies the maximum-rate user from among the Musers, and in step 1 06 
schedules data transmission to the maximum-rate user in the n-th transmission slot. The revenue 
vector is periodically adjusted as indicated in step 108 so as to compensate for observed deviations 
between actual and target throughput ratios. Step 108 may but need not be repeated for each 
scheduling operation, and is therefore shown in dashed outline in the figure. After the scheduling 
operation in step 106, and the revenue vector adjustment, if any, in step 108, the process returns to 
step 1 00 to schedule data transmission for the next transmission slot. 

In an altemative embodiment, instead of the expected feasible rate, one can take ^^(^) := 

+ RJ,n), using to obtain a weighted combination of received rates and slot allocations. More 
particularly, by choosing suitable values for X"^, one can give weight to balancing the proportion of 
slots allotted to the various users, while also achieving relative throughput targets. 

Discrete Rate Distribution 

This section will consider a case in which the feasible rates , . . . , 7?^ ) have a discrete 
distribution on a bounded set J^R^ . Since the feasible rates are assumed to be stationary in the 
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illustrative embodiment, the description will focus on the class of stationary policies. However, the 
analysis may be readily extended to deal with non-stationary policies. 

Let be the stationary probability that the feasible rate vector is j e /. Define R^^ = for j 

^ • - Jm)^^' Let be the probabiUty that policy % selects user i for transmission when the 

feasible rate vector is j E J. Then the minimum average throughput achieved under poHcy n is z"" 
= min 7^"" with = ^Pj^ij^^- • Let ti"^ be the revenue-based scheduling strategy corresponding 

M 

to the vector w = (w^, . . . , iv^). Without loss of generahty , it may be assumed that ^w^= I, since 

only the relative values of the revenue coefficients matter. 

Policy 71 is optimal if x^^ , are an optimal solution to the following linear program: 

maxz 

subz< ^ ;?^i5;x,^ i=l.,.,M 

M 
1=1 

> 0 i= l...,Mj'eJ, 

This imphes that for an optimal policy 7t at most | J^l + M- 1 of the variables are non- 
zero, which forces most of the variables to be one. Thus, only for a limited number of rate 
combinations are the slots shared among several users. 

It was noted above that a revenue-based policy which balances the throughputs is optimal. 
It can also be shown that the revenue criterion is in fact a necessary optimality condition, in the sense 
that there exists a revenue vector such that when user i does not have the maximum rate-reward 
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product, i.e., R < max w^R , then = 0 , i.e., user i should not be selected for 

m=l, ,M 

transmission. Thus, any optimal scheduling strategy must be a revenue-based policy associated with 
\v. More specifically, if poUcy % is optimal, then there exists a vector >0 such that 



for all z = 1, . . . , M, j e J, 



Continuous Rate Distribution 

This section will consider a case in which the feasible rates {R^,,..,Rj^) have a continuous 
distribution on a bounded set {7 c . 

Let p{u) be the stationary density of the feasible rate vector, i.e., the probability that the 

feasible rates are in some set Fc f/is ^p{u)du, Dofmo R^(u)=u^ for u = (uy, . . . , w^) e U, Let x^{ii) 

be the probability that policy % selects user / for transmission when the feasible rate vector isugU. 

Policy % is optimal if and only if (iff) x"^ (u) ^z"" are an optimal solution to the following 
mathematical program: 



maxz 



sub z < j p{u)R^ {u)x^ {u)du i= l,...,M 



ueU 

M 



xXu)> 0 l...,M,ueU. 
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As noted previously, a revenue-based policy which balances the throughputs is optimal. It 
can also be shown that the revenue principle is in fact a necessary optimahty criterion, in the sense 
that there exists a revenue vector w* such that if user / does not have the maximum rate-reward 
product on some set of non-zero measure, then user i should not be selected for transmission on that 
set. Thus, in the above sense, any optimal scheduling strategy must be a revenue-based policy 
associated with w\ More specifically, if policy % is optimal, then there exists a vector w* > 0 such 
that 



The previous two sections illustrate that revenue-based policies optimize throughput relative 
to pre-specified target values. However, as noted above, calculating the optimal revenue vector 
directly is a comphcated problem, requiring detailed information on the channel statistics in the form 
of the joint stationary distribution of the feasible rates {R^, . . . , The present invention avoids 
this problem by providing adaptive algorithms for determining the optimal revenue vector in an 
iterative fashion without the need for direct estimation of the frequency of occurrence of particular 
user rates. Starting from an arbitrary initial vector, the algorithms iteratively adjust the revenue 
vector coefficients to compensate for observed deviations from the target throughput ratios. The 
corrections ensure that discrepancies in throughput caimot persist. To ensure convergence to the 
optimal revenue vector, the size of the adjustments is gradually reduced. 

The next two sections will assume that the channel state is governed by a discrete time 
irreducible Markov chain with a finite discrete state space S. When the channel state iss E S, the 
feasible rates have a continuous M-dimensional distribution (■) on i? c [i?^.^^ ^maxl^? 0 < R^^^ < R^^^ 




for all / = 1, . . . , M 



Adaptive Algorithms for Adjusting Revenue Vector 
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< «5 with zero probability measure in any set of Lebesgue measure zero. In practice, the feasible 
rates will typically have to be selected from a limited set of discrete values. However, it is possible 
to adhere to the above assumptions by simply adding a small random perturbation. By choosing the 
random perturbation to be sufficiently small, the true achieved throughputs should be arbitrarily 
close to the perturbed ones. 

M 

Denote by :== w 6 { M ^ : ^ l} ^he set of all revenue vectors. For any w E W, 

denote by S^(w) the expected average throughput per slot received by user m under revenue vector 
w in stationarity. Define H.Jw) := —T^^M^E^Jw) := min S^(w), and E^^^(w) ~ 

M m=l, .,M 

max as the average, the minimum, and the maximum expected throughput per slot under 

m=\,...,M 

revenue vector w over all users, respectively. 

The above assumptions ensure that the expected throughput vector (Hi(iv), , . . , H^iv)) is 
completely determined by the revenue vector w, without the need to specify a tie breaking rule. The 
assumptions further imply that the expected throughput vector (Si(>v), . . . , H^w)) is a continuous 
function of the revenue vector w. 

Denote by w* the optimal revenue vector, i.e., the revenue vector which balances the expected 
throughputs. To simphfy the description, it will be assumed that w* is unique, although the analysis 
may be modified in a straightforward manner for the case where there is a range of optimal revenue 
vectors. 

Two-User Adaptive Algorithm 

For illustrative purposes, an adaptive algorithm will be described in this section for the case 
of two users, denoted user 1 and user 2. The next section will describe an adaptive algorithm for the 
more general case of an arbitrary number of users. 
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Before describing the algorithm for the two-user case in detail, some additional notation will 
be introduced. First, let w=-w^, so that ^2 = 1 - Denote AY(n) := Y^(n) - Y^in), and define U(N) 

N 

:= ^ A Y{n) as the difference in cumulative throughput between users 1 and 2 after N slots. The 

absolute difference \u{N)\ is referred to as the throughput gap. The throughput gap widens in the 
A^-th slot if \U{N)\ > max \u{n)\ . User 1 is said to be leading if U(N) > 0, and is referred to as 

lagging otherwise, and vice versa for user 2. A cross-over occurs in the A^-th slot if the leading and 
lagging users exchange positions, which means that the throughput gap changes sign, i.e., U(N) U(N- 
1)<0. 

The adaptive algorithm for the two-user case may now be described as follows. In every slot, 
the user with the maximum revenue-rate product, at the current price, is selected for transmission. 
Thus, the n-th slot is assigned to user 1 if w(n)R^(n) > (1 - w(n))R2(n), and to user 2 otherwise (ties 
being broken arbitrarily). To drive the revenue vector w(n) towards the optimal revenue vector w*, 
the price is adjusted over time on the basis of the observed throughput reahzations. As long as the 
throughput gap does not widen, the price is left unaltered. However, if the throughput gap does 
widen, then the price is changed in favor of the deficit user, thus at the expense of the surplus user. 
The price of the leading user is decreased by 5^(„), while the price of the lagging user is 
simultaneously increased by the same amount. 

To ensure convergence, a reset is triggered at every cross-over. The step size 5^(„) is then 
reduced by incrementing k(n), with {8^, ^= 1,2, ... } being a predetermined convergent sequence 
(e.g., 5^^i = 5/ with/? < 1, or 6^ - 5^k~^ with ^ > 1). 

For the two-user adaptive algorithm described above, the revenue vector w(n) converges to 
the optimal price w* wp 1 (where wp 1 denotes "with probabihty 1," meaning "always" or "almost 

surely") and consequently the sequence z{n) converges to the optimal value wpl. 
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Adaptive Algorithm for Arbitrary Number of Users 

The adaptive algorithm described in the previous section for the case of two users may be 
extended directly to an arbitrary number of users through modification of the rule for when to trigger 
a reset. For example, if a reset is triggered at every cross-over of any pair of users, as in the two-user 
case described above, then resets may occur too rapidly in the case of several users. Two 
leapfrogging users may cause the step size to be quickly reduced, while still far removed from the 
other users. The revenue vector may then get trapped in a bias region, and never reach the optimal 
point. A more appropriate rule for the several-user case is to trigger a reset only when every user 
has become leading or lagging. It is nonetheless important to ensure that resets occur frequently 
enough compared to wrong moves, because otherwise the revenue vector may continue to become 
trapped in a bias region. 

This section describes an alternative adaptive algorithm that is particularly well-suited for 
the case of an arbitrary number of users. This algorithm makes revenue vector updates based on 
sample periods of predetermined ever-increasing size. Thus, the revenue vector updates occur at 
predetermined slots K(n), instead of randomly determined slots as before, with L{n) K{n + 1) - 
K(n) the length of the n-th period. In every slot of the n-th sample period, the revenue vector 
sequence w{n) is used for selecting a user for transmission, with n used to index sample periods, 
rather than transmission slots as in the previous description. 

To drive the revenue vector w(n) towards the optimal point w*, the prices are adjusted over 
time on the basis of the observed throughput reahzations. The direction in which the revenue vector 
is modified at the /2-th update is determined by a random vector v(w(n)) based on the throughput 
obtained during the n-th sample period when the revenue vector w(n) is used. The size of the n-th 
update is 5(n) = 5j^„), with {5^^, ^= 1, 2, . . . } being a predetermined convergent sequence. Thus, at 
the (n + l)-th update, the revenue vector is recursively determined as 

w(n H- 1) ^ w(n) - 8(n)v(w(n)). 
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To ensure convergence, the step size 5{n) is reduced by incrementing k(n) every time a reset 
is triggered. Intuitively, resets should occur only rarely when the revenue vector is far from the 
optimal point w*, but should occur readily once the revenue vector is close to w\ 

The update direction v(w(/?)) is determined as follows. For every user, the empirical average 
5 throughput over the sample period is computed. The users are then partitioned into two groups: (a) 
those with above-average throughput; (b) those with below-average throughput. The prices of the 
above-average users are decreased, while the prices of the below-average users are increased. This 
ensures that, as the sample size grows, the empirical average throughputs will with highprobabihty 
line up with the true expected throughputs, and the revenue vector will get closer to the optimal point 
10 w\ 

i|| More formally, denote by the throughput received by user m during a particular sample 

\ =1 period in which revenue vector w is used, define X^^^ : - — ^^e average throughput over 

all users, and denote by Q' := {m:X^< X^^J and := {m : ^X^vd the groups of below-average 
and strictly above-average users, respectively; the revenue vector update direction v(w) is then 
;Ip determined as 

/ w 



20 

Note that Q~ is always non-empty, since it is impossible for all users to have strictly above- 
average throughput. However, may be empty in a case in which all users have exactly equal 
throughput. In that case, the revenue vector is simply left unaltered. 
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Also note that the price ratios within both Q' and are maintained. This ensures that the 
expected throughput of the below-average users increases, while the expected throughput of the 
above-average users decreases. 

The above revenue vector update cannot be apphed in a case in which the prices of some of 
the users in are zero. To prevent that situation from happening, the revenue vector update process 
may be restricted to the set := {we W\ > vforallm = 1, . . . , M) ,v^\thv :^ R^J{R^^^'^ {M - 
l)^max)- If t^cn H^(>v) = 0, wMch imphes that v/ e W^, In order to restrict the revenue vector 
update process to the set W^, the update is truncated at the boundary if necessary. 

To ensure convergence, a reset is triggered under the condition that every user has been a 
member of fi^ at least once during a consecutive sequence of updates. Once the reset has occurred, 
the next one is not triggered until every user has been a member ofCt at least once again. 

The above-described revenue vector update increases the throughputs of the users in £l~ and 
decreases throughputs of the users in Q"^. This can be seen as follows. Let w PFbe two revenue 

vectors and 0" 0 c {1, . . . , M\ be two groups of users such that for all / e 0" — > — for all k 
^ /, and for all / e 0', — < — for all k ^ /. Then 

As indicated above, when a correct revenue vector update direction is selected, the revenue 
vector should get closer to the optimal point w* by some definite amount. To measure distance from 
, a function r(') is introduced which attains a unique minimum at w*. Define := {>v e : 



19 



Borst 13-7 

^max(^) " ^mm(^) ^ ^} an "E-neighborhood" of w*. The following assumption implies that if a 
correct direction is chosen, then outside the reduction in the value of r(') for small enough step 
size 8 is at least 5 times a constant of proportionality r|. 

There exist positive constants 6* > 0, r| > 0, ^ > 0 such that for all revenue vectors w ^ F^, 
for any ^-right direction v(w), and for any 5 e (0, 5*), 

r(w+ 8v(w)) < T{w) ~ 5r| . 
One example of a function T(') satisfying this assumption is 

i.e., the maximum difference in expected throughout between any pair of users. By definition T(w*) 
= 0, and T(w) > 0 for all w ^ v/, with strict inequality in case the optimal revenue vector w* is 
unique. Another example of a function T(') satisfying the above assumption is 

M 
m=l 

i.e., the total expected revenue earned. As noted previously, the optimal revenue vector w* 
minimizes that quantity over all vectors in the set W, i.e., T(w* ) < T(w) for all w E W,w ^ w* , with 
strict inequaUty in case is unique. 

The second example T{') function, unlike the first, also satisfies the above assumption for 
various alternative options for selection of a revenue vector update direction, e.g., 

v* = l-P^>0 z* = arg min X^, 



20 



Borst 13-7 



V . = - 1 y = arg maxX^ 

m=l,...,M 



and = PJ{M- 2) for all k ^ ij\ for lim„ ^^ = 0,3. given positive sequence. This alternative option 
is referred to herein as the "Update-Extreme" algorithm, while the basic implementation previously 
5 described will be referred to as the "Move-to-Average" algorithm. 

The revenue vector sequence w(n) for both of these algorithms converges to the optimal 

revenue vector w wpl, and consequently the sequence z(n) converges to the optimal value z wp 
1. 

I ll The ;7-th sample is referred to as "^-right" if for every user the empirical average throughput 

;llO is within ^ from the true expected throughput, i.e., ^^(^) - S ^ (w(^))| < ^ for all m = 1 , . . , , M. 

If Otherwise the sample is referred to as "^-wrong." 

The following points can then be made regarding the above-described example adaptive 
Q algorithms: 

jS; L For any fixed ^ > 0, the total number of ^-wrong samples is finite wp 1 . 

^£5 2. The total number of resets is infinite wp 1. 

\^ 3. The revenue vector sequence w(n) cannot converge to a point outside F^. 

4. The revenue vector sequence w(n) visits F^ infinitely often. 

5. The revenue vector sequence w{n) cannot move from F^ to outside F2e infinitely often. 
The illustrative embodiment as described above establishes convergence to the optimal 

20 revenue vector w* through selection of the step sizes {5^^, A: 1, 2, . . . }as a convergent sequence. 
As one possible alternative, the step sizes may be kept fixed at a given value 5. In this case, it is 
expected that the revenue vector sequence will continue to oscillate around w*, but with smaller 
amplitudes for smaller values of 5. There is an inherent trade-off between the accuracy achieved on 
the one hand and the speed of the convergence, and thus the responsiveness to changing conditions, 

25 on the other hand. The value of 5 then may be used to find the right balance between these two 
conflicting objectives. 
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Numerical Results 

This section describes numerical results of experiments conducted to investigate the 
convergence properties of the adaptive algorithms described above. The results verify long-term 
convergence of the adaptive algorithms, and also illustrate transient performance, in particular the 
5 rate at which the revenue vector converges to the optimal vector. 

Results are presented below for a total of seven experiments. The first three experiments to 
be described consider continuous rate distributions. The fourth experiment assumes a discrete 
distribution where the feasible rates are determined by a fading process via signal-to-noise ratio . The 
fading process is modeled using a discrete number of sinusoidal oscillators as described by Jakes' 
1 0 model, W.C. Jakes, "Multipath Interference," Microwave Mobile Communications, Ed. W.C, Jakes, 
: |i IEEE Press, Piscataway, 1 974. The final three experiments examine how well the throughput ratios 
; S ; are maintained, and how well the adaptive algorithms are able to track changes in channel conditions 
Ul or throughput targets . 

,1^ The first experiment considers a scenario with two users having independent exponential 

15 rates. The feasible rate for user / is governed by a conditional exponential distribution on an interval 

20 withG; = I - e~ ' """^ a normahzation coefficient, z = 1, 2, The corresponding density /(*) for 
user / is 

25 The normalized expected throughput received by user 1 as a function of the revenue vector w is 
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and similarly for user 2, 



5 




=0 For this experiment, it was assumed that [R^^, R^J - [10,400] Kbits/s and (y^, ^ (0.02, 

13 0.01). Thus, the feasible rate for user 2 is about twice as large in distribution as for user 1. The 
iSJ throughput target for user 2 is also set twice as large as for user 1, i.e., (aj, aj) = (1, 2). 
fO FIG. 4 shows a plot of the normalized expected throughput per time slot H/w) for the above 

H parameters as a function of revenue vector w. From the figure, it can be seen that the optimal 

revenue vector is w* ^ 0.6, which may more precisely be determined as w* ^ 0.593 using bisection. 
Q In the first experiment, the two-user adaptive algorithm described previously was run for 

Pi 1000 slots, using step sizes with initial value 5j = 0.5 and reduction factorp = 0.9. It was 

^tS found that the revenue vector converges to the optimal vector in roughly 300 slots, which 

corresponds to about 0.3 seconds of operation. 

The first experiment was repeated for non-geometric step sizes 5^^ = S^yf^ withy5 successively 

chosen as 1.5, 2.0, 3.0, 4.0. It was found that the sum of the revenue vector changes is still 

convergent, although the step sizes decay slower than before. In addition, convergence is 
20 considerably slower for smaller values of y5, i.e., slower decay of the step sizes. 

The second experiment considers a scenario with three users. As before, the feasible rate for 

user / follows a conditional exponential distribution on the interval [10, 400] with parameters (yj, 
Ts)^ (0.02, 0.01, 0.02). Thus, the feasible rate for user 2 is about twice as large in distribution as 

for users 1 and 3, 
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The target throughput ratios for the three users are set equal, i.e., (a^, a3 ) = (1, 1, 1). The 
optimal revenue vector is w* - (0.424, 0.152, 0.424) as may be determined using numerical 
integration and two-dimensional bisection. It was found that the optimal prices for users 1 and 3 are 
higher than for user 1, as is required in order to obtain equal throughput, since the feasible rate for 
5 user 2 is stochastically larger. 

In the second experiment, the above-described Move-to-Average and Update-Extreme 
algorithms were run for 5000 slots, or approximately 5 seconds of operation, with L{n) = lOn slots 
for the n-th update. This amounts to roughly 30 revenue vector updates. The initial revenue vector 

is set to w(l) = (0.3, 0.6, 0.1). The step sizes were 5^ = k'^,k=l,2, The resulting revenue 

10 vector trajectories converge as expected. For example, the revenue vector for the Update-Extreme 
i algorithm after 30 updates is w(30) - (0.441, 0.123, 0.436), quite close to the optimal one. 
~{ The second experiment was repeated for the Update-Extreme algorithm using 40n and 60n 

slots for the n-th update, with the same power series for 5^. As expected, using fewer samples per 
|: update leads to a slower and 'noisier' convergence to the optimal revenue vector w\ 
rs The third experiment considers a scenario with eight users. As before, the feasible rate for 

^ user i follows a conditional exponential distribution on the interval [10, 400]. The exponents were 
□ chosen at random uniformly in [0.01, 0.05], and turned out to be approximately (0.0489, 0.0263, 
:;j 0.0139, 0.0480, 0.0220, 0.0107, 0.0461, 0.0128). 

The target throughput ratios were again set equal for all users. As before, it was expected 
20 that a larger value of the exponent y, inducing smaller feasible rates, requires a higher price in order 
to obtain equal throughput. 

The Move-to-Average and Update-Extreme algorithms were run for 15,000 slots, or 
approximately 15 seconds of operation, withZ(n) = 30n slots for the n-th update. This amounts to 
roughly 55 updates. The initial revenue vector is set at random. The step sizes were 5^ = k'^,k=l, 

25 2, The resulting revenue vector trajectories converged as expected. 

As previously noted, the fourth experiment considers a scenario with discrete rates govemed 
by independent fading processes as described by Jakes' model The mean received powers of user 
1, 2 and 3 are -15.0 dB, 0.0 dB, and -10.0 dB, respectively, which in conjunction with the Jakes' 
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model fading process led to the following feasible rates per slot as a function of signal-to-noise ratio 
(SNR) in dB: 



SNR (dB) 


Rate (bits) 


-5.0 < SNR 


1000 


-10.0 < SNR < -5.0 


500 


-20.0 < SNR < -10.0 


250 


-30.0 < SNR < -20.0 


100 


SNR < -30.0 


30 



% 

The throughput target for user 2 was set twice as large as for users 1 and 3, i.e., (a^, a^) =^ (1, 2, 
^1 1). The Move-to- Average and Update-Extreme algorithms were run for 10,000 slots, with L(n) 
%| n slots for the n-th update, using step sizes 5^ = k'^^^ and 6^ Af^, ^ ^ 1, 2, . . 

As noted previously, the discrete rate values are perturbed by adding a small uniformly 
W distributed random variable to obtain a continuous version of the problem. This ensures that the 
i|i optimal control algorithm is determined by the revenue vector only. 

: "^^f In the fourth experiment, both algorithms converged as expected. For the Mo ve-to- Average 

algorithm, the realized throughputs are reasonably close to the target ratios, provided the step size 

20 is reduced sufficiently slowly. In addition, the revenue vector converges fairly quickly, unless the 
step size is reduced so fast that the process gets essentially overdamped. For the Update-Extreme 
algorithm, the throughputs are approximately 130 bits per slot for both users 1 and 3 and 270 bits 
per slot for user 2, quite close to the target ratios. In addition, the revenue vector converges to the 
optimal value in about 5 seconds. 

25 The fifth experiment was a repeat of the previous experiment for a forcing scheme which 

assigns the n-th transmission slot to the user m\n) with the current minimum normaUzed throughput, 
i.e., 
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m (n) = arg min yjn)/a^- 

m=l, ,M 

By construction, the forcing scheme reahzes the target throughput ratios perfectly, in the sense that 
with probabihty 1, 

5 

yj{N) 

ill for all pairs of users / J = 1, . . . , M. 

: The downside of the forcing scheme is that it generally achieves lower throughput in absolute 
\W terms, as it does not take advantage of the variations in the feasible rates. Under independent 
= rj identical distribution (i.i.dO assumptions, the throughput obtahied under the forcing scheme may be 
computed in closed form. The decision as to whether or not the n-th slot is assigned to user i is 
III entirely determined by the normahzed cumulative throughputs, which only depend on the feasible 
if! rates in previous slots. Under the i.i.d. assumptions, the feasible rate for user / in the n-th slot is 
ii5 independent of the feasible rates in previous slots. 

The above-described forcing scheme was run for a period of 5000 slots. The achieved 
throughputs were approximately 90 bits per slot for both users 1 and 3, and 1 80 bits per slot for user 
2. The results show how tightly the target throughput ratios are maintained under the forcing 
scheme. In absolute terms, however, the throughput for all users is about 30% smaller than for the 
20 corresponding revenue-based scheduling technique of the present invention. 

The remaining experiments to be described examine how well the adaptive algorithms of the 
invention are able to track sudden changes in the target throughput ratios or channel conditions. Each 
of these experiments considers a scenario which includes three users. 

In the sixth experiment, the throughput target for user 3 is initially set to some low value. 
25 After 80 seconds, the throughput target is suddenly incremented to allow for the transmission of a 
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data burst for user 3. The results for this experiment show that after a few oscillations the revenue 
vector quickly settles down to the new optimal value. 

In the seventh and final experiment, the control is "cycled" approximately every 5 seconds. 
To test the tracking capability, the mean received SNR of user 3 is lowered at a rate of 5 dB/sec for 
5 5 seconds. This is expected to lead to a rapid change in w\ The change in SNR is initiated after 1 5 
seconds of simulation time, and stopped 5 seconds later. 

Both the Move-to-Average and the Update-Extreme algorithm were run, for two different 
step sizes, 5^ = k'^ and 5^ - As expected, for both algorithms the revenue vector converges more 
slowly for the former step size. More specifically, with 5^ = ^"^ convergence to the optimal revenue 
.10 vector occurs after about 25 seconds, while for 5^ = k'^^^, the optimal revenue vector is approached 
^11 shortly after 20 seconds but with stronger fluctuations around the optimal value. 

In the interval where the power is being changed, the revenue vector adjustment remains 
fairly large, which is an advantage of the reset conditions. 
J: It should be noted that the sequences L(n\ 5^, or the cycle interval used in the above- 

J5 described experiments should not be viewed as optimal values. Those skilled in the art will 
■-^ recognize that other values can be used, depending on the requirements of a given application. It 
□ should also be noted that filtering may be used to remove high-fi*equency components fi*om a control 
signal used to adjust the revenue vector. 

The foregoing description of the invention is intended to be illustrative only. For example, 
20 the particular configuration of the system 1 0 of FIG. 1 and the base station of FIG. 2 are for purposes 
of illustration only, and should not be construed as limiting the invention to any particular 
embodiment or group of embodiments. In addition, the particular data transmission scheduling 
implementations described herein are by way of example, and alternative embodiments of the 
invention may use other arrangements of processing steps. 
25 It should also be noted that the example adaptive algorithms described herein may be 

enhanced by allowing the step sizes or the sample sizes to be adapted in response to non-stationary 
changes in the feasible rate declarations. 

Since the adaptive algorithms of the present invention require only observations of feasible 
rates, they may also be used for admission control purposes. An advantage of such an arrangement 
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is that a prospective user need not be allocated any resources until the admission control decision 
has been made. 

The illustrative embodiment described above considered a scenario with only one user 
scheduled at a time and a single rate sample per user per slot. These conditions, however, are not 
essential for the underlying optunaHty principle to apply. Revenue-based policies which balance the 
throughputs continue to be optunal in situations where several users may be scheduled at a time and 
various auxiliary decisions may be taken. 

As an example, consider a throughput optimization problem for two adjacent base stations, 
one located in cell 1 and another in cell 2. Let be the rate in a given slot for user m in cell 1 if 
both base stations transmit, and let /? ; be the rate for user m if only base station 1 transmits. Let i?; 
and i? ; be defined similarly as the rate m a given slot for user / in cell 2. A revenue-based policy 
may then select the decision which maximizes revenue over all feasible options: 

w € 1, / e 2 
m e 1 
U2 

In this example, the decisions as to which users are scheduled and which base station transmits (1 
or 2 or both) are taken jointly. The revenue vector w* which balances the throughputs will be 
optimal, and may be found using the adaptive algorithms described previously. This approach may 
also be used in conjunction with scheduling of other system resources, such as antermas of an 
antexma system. 

As previously noted, the data transmission scheduling techniques described herein can be 
implemented in whole or in part in the form of one or more software programs executed using 
appropriate system processing hardware, such as processor 36 and memory 38 in the apparatus of 
FIG. 2. 



Revenue ^ 
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It should be noted that the above-described embodiments are illustrative only. For example, 
the adaptive algorithms described herein may be utiUzed for apphcations other than the specific 
communication network applications described herein. 

As one particular illustration, the adaptive algorithms may be used in off-hne applications 
5 such as cellular system planning and design. In a voice network using frequency hopped GSM vv^ith 
synchronized transmission slots, there may be, e.g., 12 frequencies and a small number of users per 
cell in each transmission slot, e.g., 3 such users, hopping between frequencies at random and whose 
signals mutually interfere. Because an encoded voice signal is made up of a number of hops, e.g., 
8 hops, the voice signal may tolerate colhsions where several users at once all randomly choose to 
10 hop onto the same frequency. The resulting interference may be controlled by appropriately setting 
fl the power levels of users. If it is desirable to achieve a particular error probability for a codeword, 
% e.g., 2 %, for the voice system, the adaptive algorithms may be used to adjust transmit power and 
find levels at which the error rate is acceptable for all users. Such an operation would typically be 
|~ performed off-line rather than in real time. 

f 5 As described previously, an adaptive algorithm in accordance with the invention can operate 

^ with multiple base stations, such as, e.g., in conjunction with a pair of base stations in a which a 
3 decision is to be made between serving one user in each base station or only one user from the two 
it base stations. As another example, consider a single cell in which K=4 transmit diversity antennas 
are being used on the downlink. Each antenna is assumed to have a distinguishing pilot and the 
20 users give some indication as to their rate if they use a certain subset of the antennas. Each user 
might therefore suggest its "best antenna" and the rate it would achieve if assigned that antenna in 
a given transmission slot. The adaptive algorithm of the invention may then be used to assign users 
to antennas, such that a pluraHty of users may be assigned per transmission slot. 

In one possible variation of the admission control aspects of the invention, an extra "dummy" 
25 version of an adaptive algorithm may be run with an extra user and its throughput target rate 
included. An admission determination may then be made by examining the rates which would be 
achieved if the extra user is included. If the result is acceptable, the extra user is admitted and the 
extra version of the algorithm is run with the corresponding new settings. Otherwise, the extra user 
is rejected and the original version of the adaptive algorithm is used. 
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It should also be noted that revenue adjustments made in accordance with the invention may 
be configured to depend on a degree of throughput difference, in place of, e.g., a constant though 
diminishing adjustment. 

Numerous other embodiments can be configured to apply the scheduhng techniques of the 
invention jointly to users at adjacent base stations, or for utilization in conjunction with diversity 
antennas. 

These and numerous other alternative embodiments within the scope of the following claims 
will be readily apparent to those skilled in the art. 
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